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<£ ; Abstract 

We consider the time delay of electron detachment from a Coulomb center and two-center systems 
in the process of ionization. It is shown that the attosecond streaking, most usual method of time 
delay measure, can be formally described by placing a virtual detector of the arrival time delay 



at a certain distance from the center of the system. This approach allows derivation of a simple 
formula for Coulomb-laser coupling that perfectly agrees with the results of numerical solution of 
the time-dependent Schrodinger equation. The dependence of the time delay upon the energy, the 
angular momentum projection, and the azimuthal quantum number is studied for the ionization 
of molecular hydrogen ion. Finally, we propose a physical interpretation of singularities, arising 
when the formal expression for the time delay is applied to the ionization of molecular hydrogen. 
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I. INTRODUCTION 

Recently the appearance of laser systems that can generate super-intense pulses as short 
as a few hundreds of attoseconds gave rise to new capabilities in studying electron dynamics 
in atoms and molecules, which were not available earlier. Among these of particular interest 
are the methods of measuring the delay of electron ejection from the atom subject to pho- 
toionization [1] . At present the measurements of ionization delay are already performed in 
noble gas atoms using the methods of attosecond streaking |2|] and RABITT (reconstruction 
of attosecond beating by interference of two-photon transition) [3j. The time delays in pho- 
toionization of molecules are a subject of growing theoretical interest (see 4| and references 
therein) . 

The numerical study [4j of the delayed electron ejection in the process of single-photon 
ionization of a hydrogen molecule H 2 revealed an interesting feature, namely, at certain 
energies of the ejected electron the angular distribution of the molecular photoionization 
time delay possesses singularities. Fig. [T^i shows the Wigner time delay tw vs the electron 
ejection angle 9 e for single ionization of oriented H 2 molecule under the condition that the 
molecular axis is parallel to the polarization of radiation, R || e. The dependence was 
calculated for different energies of the ejected electron using the external complex scaling 
(ECS) method |5(]. For the energy of ionizing radiation photon hjj = 84 eV there is a 
discontinuity at the angle 9 e = 39°, in the vicinity of which tw tend to ±oo. From Fig Jib 
one can see that this discontinuity coincides with the zeros of the ionization differential 
cross-section. Although the probability of electron ejection in the directions close to this 
zero is small, it is still possible in principle to measure the time delay for the electron ejection 
in these directions. Large negative values of tw contradict the causality principle [6] , since 
formally they mean that the electron is ejected by the molecule long before the absorption 
of the photon. 

The aim of the present paper is to clarify the conditions that cause singularities in tw 
and to resolve the paradox of large negative delay in two-center systems. We consistently 
reconsidered the physical meaning of Wigner time delay tw for Coulomb systems, studied the 
dependences of tw upon the energy and the ejection angle by the example of the ionization 
of the model two-center system (the molecular ion H^), and analyzed the dynamics of the 
ejected electron wave packet under the conditions, for which the singularities of tw take 
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FIG. 1: (Color online) (a) Wigner delay and (b) single differential cross-section -^- of the single 
ionization of molecular hydrogen H2 versus the ejection angle 8 e at R || e for different values of 
the ejected electron energy. 



place. 

The paper is organized as follows. In Section |TT] the general question of choosing the 
origin for the time delay in the case of Coulomb field. In Section |HT1 we demonstrate that the 
known method of attosecond streaking is equivalent to a device, located at a certain distance 
from the center and detecting the delay of the electron arrival, and that the Coulomb-laser 
coupling [7|, |S|, arising in the theory of the attosecond streaking, is due to the Coulomb 
advance of the electron arrival at this virtual device. In Section IIVI we study the energy 
dependence of the time delay for spheroidal Coulomb waves. Section |V] is devoted to the 
origin of singularities in the angular distribution of the Wigner time delay in two-center 



targets, and Section [VI] presents the physical interpretation of these singularities. 
Below we use the atomic units of measurement, except other declared. 

II. TIME DELAY IN ONE-CENTER SYSTEM WITH COULOMB POTENTIAL 

In the framework of quantum mechanics the behavior of a system, ionized by an external 
laser field, is described by the wave function. After the external action is finished, it can be 
presented in the form 

^(r,t) = |/(k)^- ) (r)e-^k, (1) 

where /(k) is the ionization probability amplitude, ^ (r) is the wave function of the 
continuous spectrum, describing the free particle with the momentum k at the infinitely 
large distance. 

If the potential of interaction of the particle with the center is short-range, then tp k (r — > 
oo) = (27r)~ 3//2 exp(ikr). Using the known relation 

2-7T 

lim e ikrnn ' = — \e ikr 5(n - n') - e~ ikr S(n + n')l 
r->oo ikr 

(where 8(n + n') is a delta- function), we get that for r — > oo the wave function can be 
expressed as 

VM)~/ \f(kn)\e iS +dk- \f(-kn)\e lS -dk. (2) 

Jo Jo 

Here the phases of the integrands 

k 2 
S± = ±kr + 5(±kn) 1, 

where 5(k) = arg/(k) is the phase of the ionization complex amplitude. In the limit t — > oo 
the major contribution to the integrals is introduced by the vicinity of the stationary points 
^o = ^o( r ? t), in which the derivative of the phase of the integrand is equal to zero, 

dS± 



dk 



d6(±kn) 

±r + 



dk 



- k t = 0, 
feo 



since far from the points ko the integrands are fast-oscillating. For k > only .S^ possesses 
an extremum, so that only the first term in Eq. (J2J) will contribute to the wave function, 



.k 2 

and ip(v,t) ~ f(k n)e tk ° r ~ l 2 4 . This fact can be interpreted as follows. At the point r = m 
at the moment of time 

r_ iom (3) 

f ~k + k dk [6) 

one can detect the particle having the momentum k = fcn with the maximal probability. 
Since the ratio r/k is the time, required for the arrival at the point r of the particle that left 
the center r = at the time t = and moved with the uniform velocity k, the expression 

ldS(k) dS.,, 

'" = OF = W^ (4) 

has the physical meaning of the time delay of the particle arrival at the distance r from the 
center with respect to that in the case of uniform rectilinear motion. Note that the condition 
that the particle should be ejected directly from the center is not necessary, because if the 
particle has the impact parameter b (which in the case of free motion is the minimal distance 
between the particle and the origin), then lim t ^ 00 (r/A;) = lim^oo a/6 2 + (kt) 2 /k = t. 

The interpretation of the energy derivative of the phase as the time delay was first pro- 
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posed by Eisenbud and Wigner. In 6] Winger has shown that if a particle falls from infinity, 
then after passing the center of a short-range potential the delay with respect to the motion 
in free space equals to 2ty/. In the same paper he has shown that if the potential is attrac- 
tive, then tw > only for near- resonance energy (i.e., when the particle is captured into 
a quasistationary state for a long time). For the rest values of energy (—r pot /k) < tw < 
(where r pot is the potential action radius), i.e., corresponds to the expected from the point of 
view of classical physics. For brevity below we will refer the energy derivative of the phase, 
t w , as Wigner time delay. 

One more important note should be made before proceeding further. If the external 
impact on the system is weak, so that one can calculate the ionization amplitude using the 
first-order perturbation theory, then 

/(k) = (kHz), (5) 

where \i) is the wave function of the initial state of the system, |k) = ip^ (r) is the wave 
function of the final state, describing the free electron with the momentum k at the infinity, 
w is the operator of external perturbation. 



As known, for a centrally symmetrical system the wave function is representable is the 
form of expansion 

y^V) = A7rJ2^ i5Ak) Y; m (k)Y em (r)RM(r), (6) 

tm 

where Rke(r —* oo) = -r- sm(kr — 7r£/2 + 5e) is the radial wave function of the partial spherical 
wave, 5e is the phase of the spherical wave. In this case the amplitude is expressed as 

/(k) = J2 A tm(k^ e e^Y em (k). (7) 

Em 

If the transition amplitude A^ m = An (Rk(Y^ m \w\i) differs from zero for only one £, then, 
obviously, the Wigner time delay depends only on the corresponding partial phase 

tw = ^ (8) 

In other words, if in an one-center system the transition occurs into a state with fixed angular 
momentum, then the Wigner time delay appears to be completely determined by the energy 
and the quantum number of the squared angular momentum, and is by no means related to 
the initial state and the particular form of the perturbation potential. 

If the potential is non short-range and tends to a Coulomb one, Z/r, at large distances, 
then the asymptotic form of the wave function becomes essentially different, since the phase 
in Eq.(j2]) will contain the term logarithmic in r, taking the form 

Z k 2 

S + = kr + -\n2kr + S(kn) 1. (9) 

The stationary point will be now determined by the solution of the equation 

dS + Z , , Z dS , . 

—± = r-—ln2kr + — + — -kt = 0. 10 

dk k 2 k 2 dk 

From this it is seen that the time of arrival of the electron, having the momentum k, at the 
point at the distance r from the Coulomb center with the charge Z will be equal to 

t(r) = r ---ln2kr + - + t w . (11) 

As it has been already mentioned, in systems with short-range potentials the Wigner 
time delay has simple physical meaning, namely, it is the delay of arrival of the particle, 
ejected from the force center, as compared to the case of rectilinear and uniform motion. 



To understand the physical meaning of Wigner time delay in the case of Coulomb field, the 

comparison should be made with a certain classical motion in Coulomb field. The most 

rational choice is the motion of a particle with zero angular momentum, starting from the 

center r = at the moment of time t — 0. The law of such motion can be obtained from the 

known general law of a particle motion with arbitrary angular momentum. However, it is 

simpler to use the known general formula for the law of the motion in one-dimensional case, 

r dr 
tc(r) = / --, (12) 

Jo P( r ) 

where p(r) = \jk 2 + 2Z/r is the momentum of the particle at the distance r from the center. 

Performing the integration, we get 

_ p{r)r Z p{r)/k - 1 
tc{r) - — + ¥ ln p(r)/k+l (13) 

The asymptotic form of this law is 

. r Z, 2k 2 r Z 
^(^oo) = ---ln^ + -. (14) 

Let us introduce the notion of the delay of a " 'quantum" ' particle moving in the asymptotic 
Coulomb field with respect to a classical particle with the angular momentum £ = 0, ejected 
from the center at the time moment t = 

t = ]im[t(r)-t c (r)}. (15) 

From Eqs. (TTTj) and (TT41) it is clear that the ejection delay related to Wigner time delay as 

Z , k . . 

t = t w + —\n-. (16) 

Fig. |2£i shows the Wigner time delays for the continuum of a hydrogen atom at different 
values of the angular momentum £, mutiplied by k 3 for clearness. It is seen that at small 
energies t\y — > — p- In -| for all values of £. In FigfSb the delays to are shown. At k — > 
they tend to infinity, but not as — -p-ln-|, but much slower, as 1/k 2 . Besides, t > for all 
values E and grows with increasing £. This can be interpreted as follows. The centrifugal 
potential can be considered as a short-range repulsive one, and the particle motion through 
it results in the delay with respect to the motion in the pure Coulomb field. 

Let us also consider the time delay for the particle, arriving at the detector located at the 

distance r from the center, in comparison with the case of uniform and rectilinear motion 

r Z 2k 2 r Z 

t D (r) = t + t c {r) - ^ ~ t - ^ In — + -. (17) 
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FIG. 2: Dependence of the time delay upon the energy of the ejected electron for a hydrogen atom: 
a)Wigner time delay; b) the ejection delay, Eq. (fTUj) . 

Obviously, at large r, when the logarithmic term, describing the Coulomb advance 



tcA(r) = t c {r) - 



k 



Z 



\ 2k 2 r 
In— -1 



(18) 



becomes dominant, the time delay will be negative. 



III. PHYSICAL MEANING OF COULOMB-LASER COUPLING IN ATTOSEC- 
OND STREAKING 

Apparently, the above speculations of a detector, providing direct measurement of the 
particle arrival delay, concerns a gedanken experiment rather that a real one. In real exper- 



8 



iments the delays are recorded using indirect approaches, the attosecond streaking |2[ and 
;he reconstruction of attosecond beating by interference of two-photon transition (RABITT) 



31 



. Is is interesting to relate the above discussion with these real measurements. In the 

present Section we demonstrate that the correction to the attosecond streaking time delay 

measurement, referred to as Coulomb-laser coupling [3], in fact is nothing but the Coulomb 

advance tcA{f), accumulated to the moment of essential variation of the laser field strength. 

This actually means that the attosecond streaking is exactly equivalent to a virtual detector, 

placed at a fixed distance from the center and measuring the time delay (ITT)) of the particle 

arrival. 

n 

In the method of attosecond streaking [9J the atom is simultaneously affected by the 
extreme ultraviolet (XUV) ionizing pulse with the duration of about 100 attoseconds and 
the auxiliary variable field in the form of a femtosecond IR laser pulse. The latter contributes 
to the momentum of the ejected electron, depending upon the time of ejection, which makes 
it possible to measure the ejection delay. 

Let the electron be a subject to the auxiliary field having the strength £{t). The potential 
energy of the electron in such field is 

V(r,t)=£(t)-r. (19) 

The detection of electrons is usually implemented in the direction of the polarization vector 
of the auxiliary field, since in such geometry of the experiment the auxiliary field does deflect 
the electron trajectory. Since during the action of the IR field the electron has a time to 
get far enough from the center compared to the atomic radius, the angular motion can be 
neglected. Thus, we get one-dimensional problem of an electron moving in the potential 

V(r,t) = £{t)r. (20) 

We start from the case hard to realize experimentally, but more transparent for theoretical 
consideration, the constant IR field, suddenly switched off at the time T: 

sal, = f fo - ' < T; (2D 

{ 0, t>T. 

We also assume the strength of this field £q to be very small. Note, that although the prac- 
tical detection of the momentum increment requires rather strong IR fields, no fundamental 



limitations follow from this requirement. The attosecond streaking works by virtue of the 
first-order effect surviving even for an infinitesimally weak field, so the assumption of the 
smallness of the £ magnitude is quite valid. 

At the time r the atom is affected by an ultrashort pulse of XUV radiation. As a result an 
electron is ejected with the energy E, equal to the difference of the mean energy of the XUV 
photon uq and the ionization potential /. Until the auxiliary IR field is unchanged, the energy 
will obviously be conserved, so that we can write the conservation law as E = p 2 (t)/2+£or(t). 
From the latter we can express the value of the additional momentum, acquired before the 
field is switched off: Ap = —£ r(T)/p(T). If in this formula we set r(T) ~ r^(T — r) and 
p(T) ~ p(°\T — r), i.e., instead of the exact position and momentum we use their values 
for the particle that appeared in the center at the time r and moves in the absence of the 
auxiliary IR field, then only the second-order terms in So will be changed. In the first-order 
approximation with respect to Sq we obtain 

r (o)( T _ r ) 

Both in systems with short-range potential and in Coulomb systems it is possible to assume 
p(°)(T— r) « k at sufficiently large T—r. For systems with short-range potential r^\T—T) w 
k(T — r — tw) and, correspondingly, 

Ap = -S (T-r-t w ). (23) 

Thus, in the plot of the dependence Ap(r) the time of ejection will manifest itself as a 
horizontal shift of the plot with respect to the straight line — Sq(T — r). 

In systems with Coulomb potential r(°'(T — r) « k{T — r — t£)[r(°'(T — r)]}. Since 
tn (see Eq.(fT71)) logarithmically depends upon the distance, one can use tnir^iT — r)] ~ 
td[k(T — t)]}, which yields 

Ap = -S {T -r- t D [k(T - t)]}. (24) 

Thus, the attosecond streaking works as a detector, placed at the distance r e // = k(T — r) 
from the point of free particle generation. 

We would like to emphasize the following feature of the approach leading to this conclu- 
sion. Although the particle changes its velocity in the time of the field action, starting the 
acceleration directly after the transition to the unbound state, the energy conservation law 

10 



allows expression of this change in terms of the particle behavior far from the center, so that 
the unknown details of the particle behavior near the center become unimportant. It may 
be said that until the field is switched on, the variation of the momentum is a dynamical 
variable that takes the fixed value only at the time moment when the field is switched off. 
Now consider the field strength £ (t) smoothly changing with time. The variation of the 
particle total energy in time obeys the equation ^-(t) = ^(r(t),t). From here it is easy 
to get that the first-order correction with respect to S(t) to the particle energy obeys the 
equation ^ - = %(t)r^(t — r). On the other hand, if the natural condition £(t — >■ oo) = 
holds, then AE(t — > oo) = kAp. From here we get the increment of the momentum 

A » = l *- A k- L *- (25) 

The distance r^\t — r) between the center and the particle is small at small t — r. With 
time the distance r^ unlimitedly grows, therefore, at t — r the integrand may be expected 
to contribute to the integral insignificantly. However, the subject of our interest is just the 
contribution to r^\t — r) from the shift that does not grow with time. Hence, for sure 
we assume that ^-(t — > r) — )■ 0. Below we will demonstrate that this assumption does not 
impede the applicability of this approach for realistic dependences S(t)). Since in this case 
r(°\t — t) does not contribute to the integral at small t — r, for short-range potentials we 
again can use the approximation A°\t — r) ~ k(t — t — tw)- Integrating by parts and using 
the inverse Taylor expansion, we get 

Ap = / -j-(t){t -T- t w }dt = £{r)t w - / £(t)dt 

dA 
= ~A(t) - -rr(r)t w ~ -A(r + t w ). (26) 

Here we introduced the vector potential A(t), related to the field force as S{r) = — ^jt(t), 
i.e., 

S(t)dt. (27) 

Thus we arrived at the common formula for attosecond streaking. In correspondence with 
it, the plot of the momentum of the registered electron versus time mimics the plot of the 
vector potential with the horizontal shift, equal to the Wigner time delay. 

In a system with Coulomb potential, assuming again that %(t — ?> r) — V to avoid the 
necessity to know the particle trajectory near the center, we can use the approximation 

11 



r^°'(t — t) ~ k{t — t — to[k(t — t)]}. Integrating by parts and using the fact that tr>{r 
to + tcA{r) and £c*a(0) = 0, we get 

d£ 



Ap =■-- / ^(t){t - r - t D [k{t - r)]}dt 



-A(r + t ) + ApcL- 



(28) 



Here the notation is introduced 
Ap C L = - 



00 dS 

— (t)t CA [k(t - r)]dt 



-£{T)t CA (0) 



£{t) 



dt CA [k{t-T)} 
dt 



dt 



/ S{t + T) 

JO 



1 



kdt. 



(29) 



y°)(H) 

The change of the momentum Apci, arising under thejoint action of Coulomb and laser 
field, is usually referred to as Coulomb-laser coupling |7|, |8|, |ll| . However, it is easily seen 
that if one substitutes the step function (J2"T1) for the field strength in Eq. (|2"9"|) . then Eq. (l29]) 
with the use of Eq. ( Tl2|) yields Ap C L = £otcA[k{T — r)], as a result of which Eq. (l28|) 
turns into Eq. (J24|) . That is, the correction Apci is connected with the Coulomb advance 
tcA( r eff) °f the particle arrival at the point separated by the distance r e ff = k(T — r) from 
the center, which, in turn, represents the distance, to each the particle can reache before the 
strong variation of the laser field. For a periodic field with the frequency u the characteristic 
time of the field variation has the order ~ 1/w, from which it follows that yields an estimate 
r eff ~ k/u. 

Let the radiation pulse have the form S(t) = £o(t) sin(a;t), where £o(£) is a slowly varying 
envelope function. Using the expression for the sinus of a sum, and assuming the integrand 
in Eq. (f29|) to decrease much faster, than the variation of £o(t), so that the approximation 
£q{t + t) ~ £q(t) is valid, we obtain from Eq. (l29l the following expression: 



ApcL 

Here we use the notation 



£o(t) . , \ T , s . £o(t) , v r ( v 
sinwricsO H cos(wr)i sin (a). 



UJ 



UJ 



(30) 



h\ n (a) 



where the parameter a is defined as 



cos(t) 



sin(t) 



l\ + 2a/ 1 

1 
'I + 2a/t 



- 1 



dt: 



dt, 



Zoo 



(31) 
(32) 

(33) 
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Note that Eq. (1301) may be considered as the result of splitting of the field strength £{t) 
into two terms, £(t) = £\{t) + £n(t), where the first term £\{t) = £ (r)sm[u)(t — r)] (that 
yields the first term in Eq. (l3"0l ). satisfies the condition -^(t — r) =0, which we use to avoid 
the dependence of the final result on the details of the particle behavior in the vicinity of 
the center. 

Under the condition a<l the integrals have simple expressions 



la 
L 



(a) 



r^-j 


—a 


ln--l 

a 


-7 




r~^ 


Ti 


2 


3 2 

2+ -a 2 
2 


a 


1 

~ 6 



3^ 



7 



(34) 
(35) 



where 7 = 0.57721 ... is the Euler constant. If the variation of the envelope function £o(^) 
during the period Tm = 2ir/u is small and £o(t — > 00) = 0, then the vector potential is 
approximately expressed as A(t) w -^- cos cut. With this fact taken into account, Eq. (l3"0"|) 
can be rewritten as 



dA 
Apcl = --j-tcA(r eff ) - A(r)a, 



where 



eff 



UJ 



exp 



3tt 



is the effective radius and 



Ti 



a 



-Lin (a) = —a— -a 



a 6 



2 2 
Finally, the increment of the momentum is expressed as 

Ap = -(l + a )A(r + ts) 



7 



(36) 



(37) 



(38) 



(39) 



where 



ts 



to{r. 



eff) 



1 + a 



(40) 



is the time shift. 

We recall that to(r e ff) arises from the first term in Eq.f )30p and the increment of the 
amplitude of pulse variation a appears from the second one. The equation (130]) follows from 
the decomposition of the periodic field oscillation into two components, one of them being 
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constant at the ionization time moment r and the other one changing at this moment. Sum- 
marising these facts, we can conclude, that, in full analogy with the above model example 
of a stepwise field, the time delay t D {r e jf) is recorded at the moment of the field switch- 
off, and r e ff is the distance reached by the particle at this moment. The field component, 
changing at the moment of ionization, gives rise to a distortion of this picture, increasing the 
momentum variation, synchronous with the vector potential of the external field. Against 
the background of this increased "swing" of the momentum the observed time delay ts is 
reduced. 
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FIG. 3: (Color online) Dependence of ts upon the energy of the ejected electron with the angular 
momentum £ = 1 for the ionization of a hydrogen atom (solid line) and a helium ion He + (dashed 
line). For the frequency of the IR field u = 1.55 eV. the results are calculated analytically using 
Eqs. (flu]) . (|37|) . (|38|). and (|33|) (lines); for comparison the results of solving the time-dependent 
Schrodinger equation for an atom in the field of a femtosecond IR laser pulse from |l0T | (points) 
and [3j (dotted line) are shown, as well as the result of the eikonal approximation 



From Fig. [3] it is seen that the simple analytical expression (T40l) provides complete agree- 
ment with the results [lfj and [3(, obtained from the solution of time-dependent Schrodinger 
equation for an electron affected by the ion field and the field of an IR laser pulse. 

Let us explicitly express the time shift in terms of the Wigner time 



ts 



t 



w 



1 + a(uZ/k 3 ) 



+ tlR, 



(41) 
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where the additive part of the contribution from the Coulomb-laser coupling to the observed 
time shift is 



tin 



' 2k 2 3tiluZ 

In 1 - 7 + 



(42) 



1 + a(ooZ/k 3 ) k 3 [ u ' 4 k 3 

Using the eikonal approximation, the authors of J8J also derived an analytical expression 
for tjR. Although the dominant term in this expression, logarithmically depending upon 
1/w, is analogous to the dominant term in our Eq. ( 1421) . the minor terms are not similar, 
and even the power of the momentum that enters them is different. At k — > oo the formula 
from 8] differs from our one by 2 r y/k 3 . From Fig 14] it is seen that our Eq. ( 142]) and the 
formula from [8j yield close but not identical results. For E = 10 eV the difference between 
the results amounts to 7 attoseconds, which is available for the present-time experimental 
facilities. In Fig. |3] one can see that our results are apparently closer to those obtained 
by means of the TDSE numerical solution than the result of l8j. Yet more important is 

n 

the difference in the approach and the interpretation of the results. In |8J tm was obtained 
by calculating the initial electron position r . Naturally, this makes it hard to generalize 
the results of [8] over potentials, strongly different from the Coulomb one near the center. 
We were basing on the interpretation of attosecond streaking as an equivalent of a virtual 
detector placed at the large distance r e jf from the center, which allows considering only the 
classical behavior of the electron far from the nucleus. The Coulomb-laser coupling in this 
interpretation is a record of the Coulomb advance of the particle arrival at the point r e ff. 

IV. TIME DELAY FOR SPHEROIDAL TWO-CENTER COULOMB WAVES 

In this Section and the following ones we consider the time delay of the ionization of two- 
center targets and the specific phenomena that occur in this case. As a benchmark object 
we choose the simplest molecule, the molecular hydrogen ion H;j~, for which in spheroidal 
coordinates the separation of variables in the time-independent Schrodinger equation is 
possible 

I „ £t I _i_ I y. _i_ iv I I _, iv I I „ _i_ Ix I 

Z= l - 2 '^' 2 ' e[l,oo); V = l - 2 '^' + »' 6 [-1,1]; 0G[O,2vr). 

and the solutions (Coulomb spheroidal functions) are known [15l. Ilq]. 
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FIG. 4: (Color online) Dependence of Coulomb-laser coupling contribution tin to the time shift 
upon the energy of the electron: our formula (solid line) and the formula, derived from the eikonal 
approximation [8] (dashed line). 



The two-center Coulomb continuum wave function may be presented as a sum of partial 
spheroidal waves 15 j 



^(r) = (2,) 3 %^TL (^A,0e) i l e-^T ml (^) T (^arccosr,,^ . (43) 



Here the spheroidal harmonics are introduced as 

Y lm (c, 6, ip) = S m i(c,cos9) j= — ; Y lm (0,8,(p) = Y lm (9,(p), 



'27T 



(44) 



where S m i (c, 77) is the normalized quasiangular spheroidal function 12[ , m is the quantum 
number of the angular momentum projection onto the molecular axis, I is the quasiazimuthal 
quantum number, c = kR/2 is the unitless parameter. The quasiradial Coulomb spheroidal 
function T m i(c,£,) has the asymptotic form 



(45) 



T m i(c, £ -»■ oo) = — sin I c£ + — - ln(2c^) - — + 8 lm 

where Z + is the total charge of the nuclei, Si m is the partial wave phase, entering Eq.( 

First let us consider the case, when after the ionization the electron appears in the state 
with the fixed angular spheroidal quantum numbers m and I. In analogy with the spherical 
case (i.e., Eq.flSJ)) one can introduce the Wigner time delay for a spheroidal partial wave 



t 



\Y 



dE 



(46) 
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The ejection delay to for the given Wigner time delay is calculated using Eq. (TT5|) . assuming 
Z = Z + = 2, i.e., in this case to is in fact the delay compared to the classical electron, 
ejected from the center of a helium ion He + . 

FigJSl presents the dependence of the ejection time delay upon the energy for different 
partial waves of the H^~ continuum for two values of the internuclear distance, R = 2 a.u. 
and R = 1.4 a.u. The distance R = 2 a.u. is the equilibrium internuclear distance for 
H^. The distance R = 1.4 a.u. coincides with the equilibrium internuclear distance in the 
hydrogen molecule H2, which makes it easier to compare the results for H^ and H2. 
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FIG. 5: (Color online) Dependence of the ejection time delay for partial waves of the contnuum in 
the hydrogen molecular ion upon the energy of the ejected electron : a) R = 2 a.u.; b)R = 1.4 a.u.. 

From FigfS] it is obvious that for small energies and any spheroidal quantum numbers / 
the time to tends to that of the helium atom with the same spherical azimuthal quantum 
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numbers. The cause is clear: at very large wavelength the particular form of the potential 
near the center becomes unimportant, and the main contribution comes from the region, 
where the potential is practically coincident with the one-center Coulomb one. At large 
energies essential differences from the one-center case arise. First, growing difference appears 
between the curves with the same I, but different m. For m = 1 the deviation from the 
one-center case is smaller, and to is greater, than for m = 0. Apparently, this is because 
the electrons with larger m move farther from the axis of the molecule. Moreover, the 
curve for I = l,m = quickly descends with growing E, due to which it crosses the curve 
for I = 0, and becomes lower than it. The explanation is that the spheroidal harmonic 
T 10(77,0) possesses a node in the plane, perpendicular to the molecular axis and crossing 
it in the middle between the nuclei. This node shifts the probability density closer to the 
nuclei compared with T 00(77, <fi), which has no nodes. Thus, in the case I — l,m — the 
electron moves through the regions, where the average attractive potential of the nuclei is 
larger than in the regions, where the electron with I = moves, and at large energies this 
effect dominates over the decelerating effect of the centrifugal potential (which is stronger 
at I = 1). As a result the electron leaves the molecule faster. 
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FIG. 6: (Color online) Dependence of the time delay for the patial wave with I = 1 and m = of 
the H^ continuous spectrum upon the energy of the ejected electron at large R. 

Note, that for R = 2 a.u. the effect of "advance" of the partial wave with I — l,m = 
is more pronounced. Therefore, it seems interesting to trace the variation of the curve for 
I = l,m = under the increase of R. Figure |6] shows the dependences to(E) for different 
values of R. It is seen that with the growth of R the dip in to(E) is shifted towards smaller E 
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and becomes deeper. At R > 2.6 a.u. the minimal to(E) becomes negative, i.e., the electron 
in the field of two centers, each having the charge Z = 1 (only one of them approached 
close enough) arrives faster, than in the field of one center with the charge Z = 2. A more 
consistent interpretation of this effect is the shift of the ejection point. 

Thus, the parameter t carries information about the potential in the region, through 
which the electron has travelled. Since the angular distribution of the differential cross- 
section at the infinity represents a picture of the angular distribution of these regions, the 
comparison of to for different partial waves can provide information of the target structure. 
However, in contrast to the one-center case, when the ejection of a particle with the fixed 
angular momentum £ is provided by the angular momentum selection rules (i.e., the angular 
momentum conservation law), for the spheroidal quantum numbers only the parity selection 
rule holds. Due to this fact the experimental measurement of time delays for partial waves 
is a nontrivial problem. 

V. TIME DELAY OF H+ IONIZATION 

Since in the process of single-photon ionization of a two-center target, e.g., H^~, the 
dipole transition generally brings the system to a state, containing more than one partial 
spheroidal wave, the time delay is expected to depend upon the angle. To avoid the problem 
with possible jumps of the calculated phase by 2tc (due to the phase uncertainty of a periodic 
function) and by tt (near the resonances), instead of preliminary calculation of the phase 
followed by its differentiation, it is convenient to calculate directly the phase derivative as 
the imaginary part of the logarithmic derivative of the ionization amplitude 

*-* 01) ■ w 

Figure [7] presents the angular distributions of the ejection delay to when H^~ is ionized in 
the initial ground state with the molecular axis parallel to the polarization of the ionizing 
field, for different values of energy and internuclear distance. It is seen that, similar to 
H2 (FigJTJ), the behavior of to(6 e ) demonstrates characteristic "zigzags". The plots of the 



i. 



differential cross-section demonstrate dips at the same positions |4| . At some energy {E ~ 
eV for R = 2 a.u., and E f» 40 eV for R = 1.4 a.u.) the swing of these "zigzags" of to(0 e ) 
becomes huge. 
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FIG. 7: (Color online) Dependence of the time delay upon the ejection angle for the ionization of 
molecular hydrogen ion in the ground state at e || R: a) R = 2 a.u.; b)R = 1.4 a.u.. 

Since the delays for partial waves (Fig. |5]) have no singularities, it is apparent that 
the latter arise from superposition of several partial waves. Eq. ( l4"3"j) allows the following 
expression of the ionization amplitude having angular dependence 



/(k e ) = Y, A i^- l e l5lm T lm (c,6 e ^ e ), 



(48) 



In i 



where Ai m is the amplitude of the transition to the spheroidal partial wave with the quantum 
numbers l,m. Figure [H] shows the amplitudes A[ of the dipole transition from the ground 
state of H^ to the partial waves of the continuum depending of the final energy. It is seen 
that for all energies only the partial waves with I = 1 and I = 3 are essential. 

Consider the ionization amplitude as a sum of contributions from two dominating partial 
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FIG. 8: (Color online) Dependence of the partial components of the ionization amplitude for the 
hydrogen molecular ion in the ground state at e || R: a) R = 2 a.u.; b)R = 1.4 a.u.. 



waves 



/ = A 10 e i6l ~ iir/2 T 10 + A 30 e iS3 - i3n/2 T 



30- 



(49) 



In correspondence with Eq. (l4"T|) . tw tends to infinity near the nodes f(9 e0 , eO ) = 0, if for the 
same angles -j^f(0 e o,4> e o) ¥" 0. Obviously, the sum of two functions with different position 
of nodes can be zero only in two cases: (i) the phases of partial waves satisfy the condition 
^30 — <5io = Ntc, where iV is an integer; (ii) one of the expansion coefficients for the given 
energy turns into zero, so that the sum actually consists of a single spheroidal harmonic. 
Figure [8] shows that the second of these statements is valid, i.e., there are such values of 
energy (E ~ eV at R = 2 a.u. and E = 39 eV at R = 1.4 a.u.), at which the amplitude for 
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I = 1 turns into zero. Just at these energies Fig. [7] demonstrates the singularity of to(9 e ). 

The situation, when a partial amplitude turns into zero at a certain value of energy, 
is commonly referred as Cooper's minimum, in honor of J.W. Cooper, who theoretically 



predicted this phenomenon in the photoionization of noble gases 13j). It was shown 14 1 
that in the case of H 2 photoionization the square modulus of the spherical p-wave has a 
minimum at w ~ 82 eV (energy E m 66 eV), so that in the vicinity of this energy the 
contribution of /-wave is dominant. 

In the case of H^ ionization, considered here, the contribution that turns into zero comes 
from the spheroidal wave with the quantum number 1 = 1, rather than from the spherical 
p-wave. The remaining wave is the spheroidal wave with I = 3, so it is more consistent to 
refer the phenomenon as "spheroidal Cooper's minimum" . 



A spheroidal harmonic can be expanded over the spherical ones 12] 



T lm (c, 9,0) = J2 4 m {c)Y em (9, 0), (50) 



and correspondingly, the expansion over spheroidal harmonics |48|) can be transformed into 
the expansion over the spherical ones. If c = kR/2 is not very large, the spheroidal harmonic 
T30 is rather close to the spherical one, Y30, however, T30 contains also the contribution from 
Y\o, so that in the "spheroidal Cooper's minimum" the amplitude of the spherical p-wave 
does not turn into zero. Moreover, since in the vicinity of the Cooper's minimum the phase 
difference of spheroidal waves is £30 — 5\q ~ n/2, the amplitude of p-does not turn into zero 
exactly at any energy, although its square modulus has a minimum near the energy of the 
"spheroidal Cooper's minimum". This is the difference between Cooper's minima in two- 



center systems and in one-center ones, like noble gas atoms 13] . However, since the Cooper's 
minimum for H^ with the equilibrium internuclear distance R = 2 a.u. is present at very 
small energy of the ejected electron, and T 30 (c — ¥ 0, 9, (f>) = Y 30 (9, (f>), the photoionization of 
H^ at small energies is to yield purely "octupole" electrons, i.e., those with £ = 3. 

Taking the expansion fj4*9l) for energies E in the vicinity of the energy of Cooper's mini- 
mum E c , for which by definition A 10 (E C ) = 0, and substituting it into Eq. (T4T1) . it is easy 
to obtain 

d5 m T w A' lQ (E c ) . , , , . 

^d^T^M^) 81 ^-^ (51) 

where A' 10 (E C ) = ^j^-\ E _ E ■ From this formula it is evident that the instability of the 
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dependence of Wigner time delay upon the ejection angle at the energy Ec is by no means 
related to the phase derivatives of the partial waves, but depends on the the difference of 
the phases themselves <5 30 — 5 W . As to the singularities, they arise at the ejection angles 
6 e , for which the function T 30 (c, 6 e ,<p) possesses nodes that do not coincide with the node 
9 e = 90° of the function T 10 (c, 9 e , (p). 

In order to study the angle- independent part of tw, consider tw, averaged over the 
directions of ejection 



§t w a^dtt 
tw ~ faWdtl 



El/' 

_ Ira 



lm\ 



lm ^ ' 



where a^ = 4^- is the ionization differential cross-section, and the partial amplitude fi m = 
Ai m i~ l e zSlm . From this formula it is easy to see that in the Cooper's minimum it should be 
tw = ^ff- Indeed, from Fig. [9] and particularly from Fig. [9b it is evident that t (E) near 
Ec has a maximum and coincides with the delay of the partial wave with I = 3. 

Note, however, that although t (E) is convenient for theoretical analysis, its experimental 
determination presents a difficult problem, since it requires multiple measurements with 
different orientations of the polarization vector of the "streaking" IR radiation with respect 
to the polarization of ionizing radiation. 

VI. INTERPRETATION OF SINGULARITIES IN THE ANGULAR DEPEN- 
DENCE OF TIME DELAY 

In the previous Section we determined the conditions, under which the singularities in 
the angular distribution of the time delay occur. However, their physical meaning remained 
unclear. Although the huge values of the time delay appear in the vicinity of zeros of the 
differential cross section, in principle, the time delay can be measured in this case, too. 
Large negative delays are in apparent contradiction with the causality principle, since from 
the formal classical viewpoint it means that the electron appears to be ejected before the 
external action on the molecule. 

It is common to assume that when the molecule is subject to a laser pulse with Gaus- 
sian envelop function, and the pulse is long enough, then the momentum distribution of 
the electrons, ejected as a result of ionization, is also Gaussian |4|. The Gaussian momen- 
tum distributon gives results in a Gaussian wave packet at large distances from the center. 
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FIG. 9: (Color online) Dependence of the delay to, averaged over the angles, upon the energy of the 
electron for the ionization of the hydrogen molecular ion in the ground state at e || R (black), as 
well as that for the delay of partial waves with I = 1 (red) and I = 3 (green): a) R = 2 ; b)R = 1.4 
a.u.. 

Following the Ehrenfest theorem, the behavior of the center of a wave packet (even if it is 
non-Gaussian) coincides with that of the classical particle, so that the classical interpreta- 
tion of the time delay, used, e.g., in Section UTTl is correct. But the packet center position 
itself depends on the packet shape distortions, arising due to quantum effects. Therefore it 
appears necessary to consider the problem in terms of wave packet motion. 

Let the system be a subject to a laser (XUV) pulse having the Gaussian envelope function, 
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the carrier frequency uq, and the duration T, i.e., the field strength expressed as 



£(t) = £oexp 



a i 



2T 2 



sinwo^- 



(53) 



In the case of weak ionizing field the first-order perturbation theory yields the following 
expression for the wave function, describing the wave packet of ejected electrons 

w 2 (k - k ) 2 ~ 



^(M) ~ / /On)exp 



Akr—iEt 



dk. 



(54) 



Here the initial spatial width of the packet wq = k$T, the central momentum ko = 
a/2(wo — h) (h being the ionization potential), and /(k) = (k|z|lcr s ) is the amplitude 
(J5j) of the system ionization by a continuous wave having the frequency u = k 2 /2 + 1\ in the 
first Born approximation. 

Consider the dipole matrix element (H9|) in the energy range in the vicinity of a Cooper's 
minimum. Since from Eq. floTj) it is obvious that the partial time delay determines only 
the general angle-independent time shift, while we are interested in the interpretation of the 
second term ( JBTT) . we will assume for simplicity that the derivative of the phase is zero. Near 
the Cooper's minimum A W (E) ~ A' w (Ec)(E — Ec) ~ A' 10 (Ec)(k — kc)kc, therefore, from 
Eq. (USD it follows that 



/(k ) « _( fc _ k c )kcA' l0 (E c )te l ^T 



lOV^e, <Pe. 



+ A 30 (E c )ie id ™r 



30{Ve,(Pe) 



(55) 



Apparently, at the angles, corresponding to zeros of the function T 3 o(# e ,0e), the electron 
wave packet becomes essentially non-Gaussian, and, therefore, the case of Cooper's minimum 
requires special consideration. 

Consider the electron wave function after ionization in the /c-representation 

wl(k - k ) 2 ' 



^(k) = /(k)./^exp 



W Q 



TV 



(56) 



Introduce the coefficients of its expansion 



C n = : / (p n {k - k )4>(k)dk 



(57) 



over the basis of functions 



<Pn(k) 



w 



2 n n\y/n 



H n (ak) exp 



w 2 k 2 



(58) 
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where H n {x) is the Hermite polynomial. The functions (!58|) in the coordinate representation 
have the form of expanding wave packets 



ip n {x,t) = x H n (x/w)exp 

2 n n\\/ r nw 



x 2 



2w- 



{l-it/wl)-i(5 n {t) 



(59) 



where w(t) = wo^/l + (t/wl) 2 is the packet width, growing in time, and the phase is ex- 
pressed as f3 n (t) = (n + 1/2) arctan^/w; 2 ,) — mr/2. The functions (p n (k — ko) in the coordinate 
representation equal to ipn{x, t) exp(ik r — i-^t). Here and everywhere below x = r — k t. 

If the width 1/wq of the packet momentum distribution is small and /(k) has no singu- 
larities and zeros of the order higher than the first, then only two first coefficients of the 
expansion, Cq and C\, will not be small. It is not difficult to show that in this case the shift 
of the center of mass of the packet (x) = (r) — kot with respect to the uniform motion with 
the velocity k will amount to 



,9K[ CrC (^i|a# o)] 
|Co| 2 + |Ci| 



{x) = 2 ^VM (6Q) 



Here the dipole matrix element is expressed as 



w 



(^i|x|^o) = -^e i(/3l - /3o) . (61) 



V2 
Using the asymptotic expression e*^ 1_/3 °) = 1 — iw 2 /t + 0(t~ 2 ), we get 



H^.'IW'^^yj 



T w„. (62) 



| *— '0 1 ~r [ *— ' 3_ | w | * — --0 1 ~r |^l| 

The first term in this expression grows linearly in time, i.e., actually describes the velocity 
correction 

_ d(x) V2 agcrco] 
^T - W|C | 2 + |Ci| 2 - { } 

This correction arises simply because the momentum distribution is different from Gaussian 
centered at k — fco- Indeed, Ak turns into zero if the coefficient describing the deviation 
from the Gaussian shape C\ = 0. 

The second term in Eg. ( 162]) describes the delay of the packet center arrival at the point 
r compared with the arrival time r/(/c + Ak) in the case of uniform rectilinear motion 

(x) - Akt V2w 3[QC ] 



'-WP 



k k |C | 2 H-|Ci| 2 
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(64) 



Note, that the delay twp of the packet center of mass arises due the temporal variation of 
the packet shape, i.e., actually, because of the fact that the common phases of the wave 
packet without a node, (3 (t), and the packet with a node, fli(t), have different time depen- 
dence. The functions (15T?1) are solutions of the time- dependent Schrodinger equation with 



the Hamiltonian of a harmonic oscillator with variable frequency w{t) = l/yw({), obtained 
from the equation for a free particle by transformation to the time-scaled coordinate sys- 



tem [l8| x/w(t). The phases (3 n (t) = —mr/2 + J Q e n (t)dt are integrals of the eigenenergies 



10 

e n = (n + l/2)w of harmonic oscillator with variable frequency 1/ \Jw{t) over time. There- 
fore, we can conclude that the delay twp is a consequence of the specific quantization that 
arises when the free particle motion is considered in the time-scaled coordinate system. 

Consider the case, when the central frequency of the ionizing laser pulse coincides with 
a Cooper's minimum, i.e., when k = he- Then 

C = -A 30 (E c )ie iS ™r 30 ; (65) 

Ci = -jL-k c A' 10 (E c )ie^r w . (66) 



The corresponding combinations of coefficients that enter the expression ( 1631) and (1641) are 



|C | 2 + |d| 2 = \A 30 (E c )r 30 \ 2 + ^_\A' 10 (E C )T 10 \ 2 ; 

C?C = -^-A 3 o(E c )A' w (E c )e l ^°- s ^T 1 oT 3 o. 

Eq. ( 164]) at w Q — > oo yields the expression that does not depend upon w and completely 
coincides with the second term in Eq. (15 ip . The shift of the central momentum at Wq — > oo 

A k = — — — — — - cos {5 30 - 5 W ), 67 

Wo T so A 30 (E C ) 

i.e., for wo = oo the velocity shift turns into zero. Near the node of the angular distribution, 
i.e., at T 30 — > 0, the wave packet center time delay 



2w 2 A 30 (E c ) . ( . . ,T 



twP = -pr a, en \ sin (^ 



30 



30 — "10 ) 



k 2 c A' l0 (E c ) VdU 1U/ T 10 - 

Thus, the time delay in the vicinity of the node of the cross-section does not tend to infinity 
in reality. On the contrary, for any final wq its value tends to zero, or, more rigorously, to 
the delay of the partial spheroidal wave with 1 = 1, neglected in the above consideration for 
simplicity. 
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The maxima of \twp\ are attained at \Cq\ = \C\\. The maximal possible delay or advance 

is 

w sm(S 30 - Sio) sm(5 30 -5 w ) 

m nf ' wl = Y — 7! — = — 71 — ' (68) 

i.e., it apparently does not exceed the uncertainty of the arrival time of the ionizing laser 
pulse, equal to its duration T. Thus, there is no contradiction with the causality principle. 

VII. CONCLUSION 

It is demonstrated that the known method of attosecond streaking is equivalent to a 
device, placed at a certain distance from the center and recording the delay of the electron 
arrival. We show that the laser-Coulomb coupling, arising in the theory of attosecond 
streaking, is related to the Coulomb advance of the particle arrival at this virtual detector. 
The energy dependence of the time delay is studied for partial two-center spheroidal Coulomb 
waves. It is shown that for the energies, coincident with those of the Cooper's minimum 
(in which the amplitude of one of the partial waves of the ejected electron turns into zero), 
at angles, coinciding with the nodes of the angular spheroidal function having the quantum 
number / = 3, the angular distribution of the Wigner time delay for two-center targets has 
singularities, in the vicinity of which the time delay takes large positive and negative values. 
By analyzing the dynamics of the wave packet of the ejected electron it is shown that these 
singularities do not lead to violation of the causality principle, since in the process of ionizing 
the molecule by an external laser pulse with finite duration the real delay or advance does 
not exceed the duration of the laser pulse. 
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